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Vou. 30. No. 248. Aveust 1944. 


THE APPLICATION OF FREE ENERGY EQUA- 
TIONS TO THE STUDY OF THE SYNTHESIS 
OF HYDROCARBONS FROM CARBON MON- 
OXIDE AND HYDROGEN.* 


By W. W. Myppzerton, D.Sc., F.R.L.C., F.Inst.Pet. 


Synops3. 


Equations have been derived relating free energy of reaction to tempera- 
ture in the synthesis of a series of gaseous hydrocarbons from carbon mon- 
oxide and hydrogen. 

In building up these equations molecular heat-tem ure curves were 
drawn for each reactant and product and molec heat-temperature 
equations were fitted over the range 25-250° C., existing equations being 
generally unsatisfactory for the temperature range covering synthesis. Some 
use has been made of molecular heat and entropy values calculated from 
statistical and spectroscopic data. 

The suppression of reactions in which carbon dioxide is formed, in spite 
of free energy changes favourable to them, is attributed to special conditions 
prevailing at the catalyst surtace. 

The oes of traces of ethylene amongst the products of synthesis 
instead of appreciable quantities predicted by the free relationshi 
is accounted for by the incorporation of ethylene actually formed into 
building up of hydrocarbon chains. 

The possibility of increasing the yields of isobutane and isobutene is con- 
sidered worthy of attention. The formation of aromatic hydrocarbons 
requires a special orientation of reactants on the surface of the catalyst, and 
the probability of the necessary conditions occurring is slight. 


Tue Gaseous HYDROCARBONS. 


The symbols employed are in accordance with the Joint Report of the 
Chemical, Faraday, and Physical Societies, 1937. 


Tue Free ENERGY EQUATIONS AND THEIR USEs. 


In a study of the thermal characteristics of a reaction such as nCO + 
(2n + 1)H, = C,Hon,., + H,O the thermal data available may be 
assembled in the form 


AH, = AHp + AT + BT? +CT?. . . . (1) 


where AH, is the change in heat content during reaction at T° K. and 
AH» may be regarded as the change in heat content in hypothetical 
reaction at absolute zero. A, B and C are constants for a reaction 
producing any one hydrocarbon. ° 

The application of equation (1) gives information necessary for the 
provision of adequate means of controlling the temperature in an exo- 
thermic reaction. Its general usefulness can be extended by deriving 
from it an equation showing the change in free energy during the reaction. 


* Received 17th February, 1942. Released by Censor May 1944. 
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We thus obtain a thermodynamical approach to a series of practically 
important questions such as the following :— 


1. What operating conditions can most profitably be explored / 

2. Is any other product or any other type of reaction of the same 
order of thermodynamic probability as the main reaction we are 
considering? When this is the case there is always a danger, to be 
explored and provided for, of some alteration in working conditions 
or in the form of catalyst suddenly provoking a disturbance and 
perhaps producing undesirable products. 

3. Can the non-occurrence of equally probable reactions be accounted 
for by a different mechanism at the catalyst surface, such that the 
conditions necessary for its accomplishment cannot be fulfilled in ‘he 
prevailing circumstances? The answer to this question may serve 
to evaluate the risk of undesirable complications dealt with in the 
previous question. 

4. Do thermodynamical considerations show that some of the more 
desirable products can be formed in larger proportion ? 


Such questions as these will be discussed with reference to the gaseous 
hydrocarbons which normally accompany liquid and solid hydrocarbons 
during synthesis from carbon monoxide and hydrogen. 

For the development of equation (1) the constants A, B and C are 
derived from molecular heat-temperature equations for all the substances 
involved in the reaction. For each reactant and product an equation is 
obtained in the form C, = a + bT + cT*. Then by algebraically sum- 
ming the molecular heats of the products and subtracting the sum of the 
molecular heats of the reactants we obtain AC, for the reaction. 


Substituting in the equation oan = AC, and integrating we obtain 
ing oT / ing 


A, B and C in equation (1). The integration constant AH» is evaluated 

by finding the value of AH, for some fixed temperature in the range for 

which the molecular heat equations are valid. This is readily done by 

employing the data on heats of combustion of the substances involved in 

the reaction. 

. Substitution from equation (1) in the thermodynamical equation 
AG 


=e =— 1 followed by integration gives us the final free energy 
equation 


AG? = + A'TinT + BT? + CT? +...+4+I17 . (2) 


The integration constant I can be evaluated in relatively few ways. 
The value of AG; must be found for some fixed temperature in the specified 
range. For this purpose the equation AG? = — RT In K is not generally 
applicable to the reactions we are considering because the equilibrium 
constant K is not often available. The most useful procedure is to apply 
the relation AG; = AH} — TAS} at T = 298-1° K., where AS> is the 
sum of the molecular entropies of the products at 298-1° K. less the sum 
of the molecular entropies of the reactants. 
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Previous thermodynamical treatment of the subject } has been marred 
by one or more of the following factors :— 


A. Inaccurate thermal data such as heat of combustion, molecular 
heat at constant pressure, and entropy. 

B. Lack of information on rates of reaction. 

C. Lack of information on the difference between free energy 
change as normally calculated for conditions in which each reactant 
and each product is at atmospheric pressure, and, under the actual 
conditions of synthesis, when partial pressures are far removed from 
atmospheric. 


These difficulties have been considered in the order given. 


A.l. The Heat of Combustion of Reactants and Hydrocarbon Products. 
In recent years highly accurate data on heats of combustion have been 


’ provided by Rossini, Jessup and their co-workers. 


A.2. Molecular Heat at Constant Pressure. 


We are concerned with the conditions in the catalyst space, and mole- 
cular heat values should therefore refer to gases and vapours. When we 
examine recorded experimental figures, we find very serious discrepancies, 
even in the case of such gases as the oxides of carbon. A critical study of 
the methods employed in the observations reveals that, where possible 
sources of error have been reduced and when special refinements in 
apparatus or technique have been introduced, the values obtained are in 
good agreement with those calculated from the characteristic vibration 
frequencies of the molecules by the application of the Planck—Einstein 
equation.” 

For such gases molecular heats have been calculated from what are 
considered to be the most reliable vibration frequencies recorded in the 
literature. 

Calculation is as yet unsatisfactory in the case of the more complex 
gases. Not all of the characteristic frequencies of a complex molecule 
can be derived from the infra-red and Raman spectra, the usual source 
of such information, and, in fact, experimental specific heat values have 
been employed to predict missing frequencies in ethylene * and ethane.‘ 
Another kind of difficulty is found in ethane and in general in molecules 
containing a C-C bond. There is undoubted evidence of a restriction in 
the freedom of rotation about this bond, and until its magnitude has 
been ascertained we cannot calculate the molecular heat due to the internal 
rotations of the molecule.*. ® 

For hydrocarbon gases other than methane, therefore, the best experi- 
mental values have been selected from the literature. 


Method of Calculation. 
The contribution made to the molecular heat at constant volume by . 
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each mode of vibration of the molecule was found by applying the Planck- 
Einstein equation in the form 


- 

X (Degrees of freedom) 
1) 
where R is the molecular gas constant, taken here as 1-9863 cal.; e is the 
base of natural logarithms; -0 is the characteristic temperature of the 
vibration = 1-4327 x (Number of wave-lengths per centimetre). 

The selection of characteristic frequencies from the spectra and the 
assignment to them of their appropriate degrees of freedom has become 
a highly specialized branch of science. An extensive examination of the 
data recorded in the literature has been made. 

The total molecular heat at constant volume was obtained by adding 
the effects due to translation, rotation and vibration. 


Type of molecule. C, = “pans. + Cros. + 
Diatomic. C,= 3R+R+ 
Triatomic linear. 

co, type. C, = +R + 
Triatomic non-linear. 


H,0 type. C, = + + Cy. 
Polyatomic. CH, type. C, = + 3R + Cy. 


Calculation of C, at atmospheric pressure. 


Conversion of C, to C, at infinite dilution was obtained from the 
equation C, = C, + R. 

The correction for increase in pressure to 1 atmosphere was made in 
the usual manner by assuming that the gas obeys Berthelot’s equation of 
state. Then if P, and 7’, are the critical pressure and temperature « 


A.3. Molecular Heat Equations. 


Many of the molecular heat equations in common use are constructed 
for temperatures beyond the range in which we are interested. New 
equations have therefore been obtained fitting the graphs molecular heat- 
absolute temperature for each gas over the range 25° to 250° C. 


A.A. Evaluation of the Integration Constant AH 5. 

The value of the integration constant AH, in the free energy equation (2) 
can be found from equation (1) when we find AH, at some definite tem- 
perature within the specified range for which the molecular heat equations 
have been constructed. 

Heats of combustion afford the means of doing this, and as they are 
usually reported for 298-1° K., this temperature is used for determining 
AH,. Reported values account for any water formed in the liquid state. 
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We must therefore take into account the latent heat of vaporization of 
water at 298-1° K. 


A.5. Molecular Entropy. 


The experimental determination of entropy involves the measurement 
of specific heat down to temperatures preferably below 10° K., and there- 
after a graphical extrapolation to absolute zero. The method of extra- 
polation has been simplified by Kelley, Parks and Huffman.* The latent 
heat appropriate to any change in state or form between absolute zero 
and 298-1° K. must be measured and used in the calculation. 

The difficult procedures in the determination make it desirable to com- 
pare experimental values with those calculated for the simpler types of 
gas from statistical and spectroscopic data. 

Important discrepancies between the two sets of values in the cases of 
carbon monoxide, water vapour and hydrogen require investigation before 
we decide to adopt any particular figure. 


Entropy of Carbon Monoxide. 

If the assumption is made that the molecule behaves as though it were 
unsymmetrical in the crystalline solid, the entropy due to rotation is 
11-375 entropy units, and the total entropy of the gas at 298-1° K. is 
47-298 E.U. The experimental value reported by Clayton and Giauque ? 
is 46-22 E.U. 

There is some independent experimental evidence that the molecule 
behaves as though it were symmetrical. If this were the case, the 
rotational entropy would be 9-998 E.U. and the total entropy at 298-1° K. 
would be 45-92 E.U. The value taken for the calculation of free energy 
is the mean value of 46-07 E.U. 


Entropy of Water Vapour. 

The calculated value is Sj, = 45-09 E.U. Direct experimental 
methods according to Giauque and Stout ® gave Sj, = 44:28 E.U. A 
suspicion of randomness in the ice crystals led to an indirect method in 
which Giauque and Archibald ® studied the reaction Mg(OH)? = MgO + 
H,0. This method gave 83,,, = 45-10 E.U. 


Entropy of Hydrogen. 

The experimental determination by Giauque 7° gave S},,., = 29-69 E.U. 
From statistical and spectroscopic data Giauque '! calculated S},,, = 
31-23 E.U. The discrepancy has been accounted for by Moelwyn Hughes 1” 
as being due to a difference between the entropies of ortho- and para- 
hydrogen. The entropy of the equilibrium mixture is calculated to be 
29-59 E.U. from statistical and spectroscopic data. 

The evidence at first sight appears to favour the adoption of the lower 
figure. It is, however, improbable that the same equilibrium mixture is 
concerned in synthesis. Furthermore, we can derive a value for the 
entropy of hydrogen from the free energy of the reaction H, + 40, = H,O 
for which four different methods give AG, = — 54,507 cal. as an 
average.13 Assuming for water vapour S>,,, = 45-10 E.U. and for oxygen 
Soop, = 49-07 E.U. we derive for hydrogen S},,, = 31-79 E.U. 
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The conditions under which the free energy of formation of water was 
determined appear to exclude the possibility of normal equilibrium mixture 
of ortho- and para-hydrogen being involved. 

The value 31-23 E.U. has been adopted. 

The entropies of gases are usually referred to standard conditions at 
298-1° K. and 1 atmosphere pressure and, furthermore, the gas is assumed 
to be in the ideal state in which the fugacity is also 1 atmosphere. When 
the gas is not ideal a correction can be applied by assuming, for example, 
Berthelot’s equation of state. A gas obeying this law has an entropy 

27P 
lower than that of an ideal gas 14 by R 30P, ( , 

The correction is negligibly small for the gases and pressure and tem- 
perature conditions we are investigating. 

For hydrocarbons other than methane the best experimental erttropies 
have been adopted. 


B. Rate of Reaction. 

A method has been described for the measurement of the time of con- 
tact, and hence the rate of reaction under the conditions of gas flow of 
practical synthesis.15 Observations will be given on the rate of formation 
of methane and carbon dioxide. 


C. Correction to be Applied to Free Energy of Reaction for Differences in 


Pressure. 


Free energy equations such as (2) are based upon the assumption that 
each gaseous product and reactant exerts a pressure of 1 atmosphere. 
We can estimate the approximate correction to be made for the actual 
conditions operating during synthesis. 

For this purpose it has been assumed that blue water gas is used for 
synthesis and that 75 per cent. of the available carbon monoxide is con- 
verted into a single hydrocarbon at 200° C. and at a total pressure of | 
atmosphere. The partial pressures of reactants and products in the 
residual gas have been calculated for these conditions. 

From the equation 


where the suffix 1 relates to the standard pressure of 1 atmosphere and x 
to the actual condition, we obtain 


P, 
G, — G, = RT =} 
and for the general reaction 
nCO + (2, + Hg, +2 nH,0 « 
we have 
G, — G, = RT in an + 2) Xx p”"(H,0) (6) 


p*(CO) x **(H,) 
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For the conversion of 1 gram molecule of carbon monoxide in a reaction 
of this type the following values of G, — G, have been calculated :— 


Hydrocarbon formed. Correction, calories. 
Ethylene ‘ 2005 
Propylene. ‘ 2092 
Butenes 2127 


There is an alternative type of reaction in which carbon dioxide is 
formed instead of water. 


2nCO + (n + 1)H, . . (6) 


For this type of reaction the corrections are :— 


Hydrocarbon formed. Correction, calories. 
Methane ‘ ° 1413 
Ethane . ‘ 1093 
1090 
Butenes ‘ 1010 


The gases-have been assumed to obey the ideal gas laws, but the 
replacement of partial pressures by fugacities makes no significant altera- 
tion at the low pressures concerned. 


Tables of Data Employed. 


Taste I. 
Statistical and Spectroscopic Data. 


Gas. 6. Reference. 
Carbon monoxide 1 
Hydrogen 
Carbon dioxide { 
1 
Water vapour . : 
3 
Methane . 
3 
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Taste II. 
Molecular Heat at Constant Pressure, 1 atm. 
Constants for equation C, = a + bT + cT® + dT®. 


Refs. for 
Gas. a. 6 x 10%. | ¢ x 10%. | d x 10°. | expt. values 
used. 
Carbon monoxide 6-936 — 001 — 0-08 + 1-625 
Hydrogen * 6-956 + 0-07 
Carbon dioxide . ° 5-79 +11-0 — 20 — — 
Water vapour . 9197 — 53 + 636 +3-0 
Ethane. 3-80 +30-0 — 21 
Propane . 8-512 +31-375 oe 22 
Butane (n) ‘ . 10-89 +40-0 + 20 — 23 
Butane (iso) 11-95 +30-0 +20-0 23 
Ethylene . 3-9 +22-0 — 21 
Propylene ‘ 3-972 +37-23 23 
Butene (1) ‘ : 4-61 +513 | — — | 23 
Butene (iso) assumed . 4-61 +51-3 — 


* In calculating the molecular heat of hydrogen the correcting term 0-000,08T was 
added.** 
+ Wilson ** has shown that correction should be made here for distortion of the 
molecule by adding 2pRT. For water p = 2-04 x 10-5 and for methane p = 1-72 x 
*10°5. The corrections have been applied. 


Taste III. 
Standard Molecular Entropy at 298-1° K. 
_ 
Gas. obs. Refs. calc. Refs. adopted. 
Carbon monoxide | 46-22 E.U. 7 45-92 E.U. — | 46-07 E.U. 
47-30 * 
Hydrogen . . 29-69 10 31-23 ll | 31-23 
29-59 ¢ 12 
Carbon dioxide . 51-11 26 50-83 — 51-00 
51-07 t 27 
Water vapour. 45-10 9 45-09 _— 45-10 
44-28 8 
Methane § . ° 43-95 28 44-46 31 44-48 
43-39 29 44-64 32 
44-60 30 44-35 33 
44-34 
Ethane 54-85 34 — 54-85 
Propane. 64-70 35 — 64-70 
Butane (n) . : 75-80 { 36 75-80 
Butane (iso) 70-00 } 70-00 
Ethylene . - | 52-48 37 | 52-48 
Propylene . ‘ 63-10 38 — -- 63-10 
Butene (1) . ‘ 72-50 { 39 } 72-50 
Butene (iso) 69-00 69-00 


* Assuming molecule asymmetric in solid. 

+ Corrected for equilibrium mixture ortho and para forms. 

} Different moment of inertia taken. More accurate method of calculation 
employed. 

§ The first two experimental values quoted contained integration of Eastman’s 
equation for molecular heat of methane. This equation is not accurate for the 
temperature range. Correction of second value gives the fourth value. 
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Taste IV. 
Heat of Combustion Sughyet, 298-1° K., 1 atmos. 
a Heat of combustion, 

: Gas. g. cal. per g. mol Refs. 

Carbon monoxide 67,636 40 

Hydrogen (H,O liq.) . 68,318 41 
Hydrogen (H,O vap.) 7 ‘ 57,877 * 42 
Methane (H,O liq.) . 212,790 43 
Ethane (H,0 liq.) . 372,810 44 
Propane (H,0 liq.) 530,570 45 
Butane (n) (H,0 liq.) 687,940 44 
Butane (iso) (H,O liq.) . . 686,310 45 
Ethylene (H,0 liq.) . 337,280 46 
Propylene (H,O 491/820 46 
Butene (1) (H, O liq.) 649,660 46 
Butene (iso) (i 


hq. 647,200 41 


* From latent heat of vaporization. 
The latent heat of vaporization of water at 298-1° K. has been taken as 10,44] 


cal. per 
TaBLe V. 
the Free Energy mee for Reactions of the Type 
2 x nCO + (2n + 1)H, = C,H,,,, + nH,O (25-250° C.). 


General equation AG; = AHo + +ATinT + BT? + CT? + DT7...+ 17. 


Hydrocarbon | B x 10%./C x 10%./D x 10".) I. 
Methane | — 45-560} 14-107 | — 3-90 | —108 | —12 | —41-680 
Ethane |— 76-249/ 26-458 | — 9-54 | -—2-15 | | —67-690 
-— Propane .|—10-13 | 33-397 | — 7-51 | —3-22 —3-4 | —61-486 
U. Butane (n) —144-05 | 42-671 | — 911 | —463 | —46 | —69-913 
Butane (iso) —145-708 | 41-610 | — 4-11 | —7-63 | -—46 | —59-081 
Ethylene. || — 45: 212 | 19-402 | — 5:57} -—215 | —23 | —51-643 
Propylene . | — 81-547} 30-981 | —10-47 | ~3-22 | —3-4 | —73-165 
Butene (1) . | —114-702 | 41-994 | —14-79 | —4:30 | -—46 | —90-010 
Butene (is0) . | —117-162 | 41-994 | —14-79 | —4-30 | —46 | —86-510 

Tasie VI. 
Free Energy Equations for Reactions of the Type 
2nCO + (n + 1)H, = + nCO, (25-250° C.). 
Hydrocarbo AH 

| A. B x 10°. | C x 108. |D x 10". I. 
Methane . . | — 54-955| 17-494 | —12-09 | +0-32 | + 2-7 | — 50-950 
Butane (n) . | —181-651 | 56-218 | —41-87 | +0-90 | +11-0 | —106-89 


Relative Values of Terms Involving B, C and D. 
on For the equations in Tables V and VI terms in D contribute less than 
: 0-1 per cent. to the total free energy. D and C terms together contribute 
he —— 2 per cent. D, B and C together contribute less than 10 per 

cent. 
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Tasie VII. 


Molecular Free Energy of Formation of Gaseous Hydrocarbons from Carbon 
Monoxide and Hydrogen. 


A. Reaction in which water vapour is formed. 
B. Reaction in which carbon dioxide is formed. 


| 


Hydrocarbon. | Reaction. | g-cal. | gical. |  g.-cal, 
Methane . A | — 34,401 | —25,165 | —25,165 
Methane . B | — 41495 | —30,738 | —15,369 
Ethane A — 52,395 | —33,540 | —16,770 
Propane . ‘ . et A — 72,493 | —43,938 | —14,646 
Butane (n) ~| A — 93,355 | —55,338 —13,834 
Butane (n) B —121,730 | —77,682 | — 9,710 
Butane (iso) A — 93,254 —54,373 | —13,593 
Ethylene ‘ A — 28,209 — 14,595 — 7,298 
Propylene ‘ , of A — 51,756 | —28,582 | — 9,527 
Butene (1) . ‘ wil A — 71,641 | —38,699 | — 9,675 
Butene (iso) | A — 73,059 — 39,505 | — 9,876 


* Per g.mol. CO converted. 


Discussion of Free Energy of Reaction. 


The high negative values of the molecular-free energy of the reactions 
and the slope of the free energy vs temperature curve show that the lower 
limit of temperature will be determined by the rate of reaction. A more 
active form of catalyst would allow of lower working temperatures. 


Formation of Carbon Dioxide. 


The reaction in which carbon dioxide is formed has a higher negative 
free energy than the reaction producing water. It would therefore be 
expected to accompany the reaction in which water is formed. In the 
normal working of the synthetic process with nickel and cobalt catalysts 
the production of carbon dioxide is suppressed, but it can be initiated 
and allowed to replace the water-forming reaction by raising the tem- 
perature five to ten degrees above the optimum point. In the intermediate 
stages both carbon dioxide and water are produced, and methane appears 
in amount about equivalent to the carbon dioxide formed. With the 
entire suppression of water formation, methane is practically the only 
hydrocarbon to be made. 

The elimination of the carbon dioxide reaction cannot be accounted for 
by a low rate of reaction. Aicher, Myddleton and Walker 15 have shown 
that on a nickel catalyst at about 180° C. the rate of conversion of carbon 
monoxide, present in water gas, into higher hydrocarbons and water is in 
the region of 20 normal litres per litre of catalyst space per hour. 

A quantity of the same catalyst in an aluminium tube controlled at 
180° C. by an electrically heated and controlled aluminium block produced 
methane and carbon dioxide from the time blue water-gas was admitted 
for an observed period of ten days. The flow of gas was maintained at 
the rate of approximately 100 litres per hour for a catalyst volume of 
200 c.c. The inlet and outlet gases averaged :— 
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From 100 vols. inlet. 
Inlet, Outlet, 
.0/ 
vol.-%. Consumed, | Produced, 
vols. vols. 
CH, 2-4 36-7 on: 19-5 
H, . . 48-4 15-4 39-2 van 
N, (uiff.) 5-5 9-9 
| 100-0 100-0 


The contraction in volume was 40-3 per cent. 

One litre of catalyst space thus converted 200 litres of carbon monoxide 
per hour into methane and carbon dioxide. 

The catalyst was found to contain about 0-5 per cent. free carbon after 
the reaction. 

These facts point to different mechanisms in the two types of reaction. 
A similar conclusion was reached by Craxford 47 from a study of the rates 
of formation and reduction of cobalt carbide, assumed to be formed as an 
intermediate during synthesis with cobalt catalysts, by reaction of the 
metal with carbon monoxide and hydrogen. 


¢ H, 
Co—Co—Co —> Co—Co—Co + H,O 


The formation of higher hydrocarbons and water is accounted for by 
reaction between the carbide and molecular hydrogen. The formation of 
methane is assumed to be due to reaction of the carbide with chemi- 
sorbed hydrogen, while carbon dioxide is produced by the process CO + 
H,O = CO, + H,. When higher hydrocarbons and water are formed, 
long-chain molecules are built up from methylene groups, the immediate 
reduction products of the carbide, and are broken down into lower- 
molecular-weight hydrocarbons by chemi-sorbed hydrogen. 

On this basis the behaviour ofNiron catalysts is anomalous in many 
ways. When they form higher hydrocarbons they produce carbon dioxide 
simultaneously to the exclusion of water.** Furthermore, the hydro- 
carbons have a much higher average molecular weight than those formed 
on nickel and cobalt catalysts. 

It appears possible to account for the behaviour of all three catalysts 
by assuming that when carbon monoxide is adsorbed by dual attachment 
at adjacent atoms on the catalytic surface, water is formed by subsequent 
reaction with molecular hydrogen. Methylene groups are formed at the 
same time, and give rise to hydrocarbon chains. When carbon monoxide 
is adsorbed by single attachment, molecular hydrogen produces methylene 
groups and carbon dioxide. The methylene groups may be converted 
into methane by the action of chemi-sorbed hydrogen. 

The difference in spacing of the metallic atoms in iron catalysts may 
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account for absence of dual attachment of carbon monoxide in this case. 
Experimental evidence is being sought. 

The view has been expressed by H. 8. Taylor ® that at 270° C. carbon 
monoxide attaches itself to nickel catalysts in the form of Ni,C and NiO. 
For the purpose of the present discussion it is not necessary to assume 
that the bond energy between carbon and nickel is sufficiently high to 
give rise to carbide. The carbon would be more reactive if this were not 
the case. 


Formation of Ethylene. 


On purely thermodynamical grounds we should expect much more than 
the traces of ethylene normally produced. Typical figures obtained during 
the treatment of 3000 cu. ft. of blue water-gas per hour at approximately 
200° C. are :— 

Carfon converted to ethylene — — — trace 
Ratio Casten to propylene 13 
Carbon to propane 
Carbon to butene 


Ratio Carbon to butane 


The free energy values point to a decline in the ratio from C, to C, with 
an appreciable formation of ethylene. The difference in the rates of 
hydrogenation of propylene and ethylene is insufficient to explain the 
absence of ethylene. We therefore seek the explanation at the surface of 
the catalyst. 

It was demonstrated by Smith, Hawk and Gold®n © that when ethylene 
is admitted along with carbon monoxide and hydrogen, a substantial part 
of it is used up in the formation of higher hydrocarbons. Twigg and 
Rideal 5! showed that ethylene is adsorbed on the surface of a nickel 
catalyst by dual attachment on adjacent nickel atoms in the 110 and 111 
planes of the face-centred cubic crystals. Herington therefore accounts 
for the incorporation of ethylene in the synthesis of higher hydrocarbons 
by assuming that synthesis normally occurs on the same planes.™ 


Formation of Branched-chain Hydrocarbons. 

The free energy relationships point to the formation of appreciable 
quantities of these hydrocarbons.. Schaarschmidt and Marder ® found 
33 per cent. iso-paraffins in the paraffinic constituents of the synthetic 
gasoline fraction, and Koch and Hilberath * reported analyses showing 
presence of such compounds as 2-methyl butane, 2-methyl pentane, 
3-methyl hexane, 3-methyl butene-1, 3-methyl pentene-1, 4-methyl 
pentene-1, 4-methyl hexene-1, and 5-methyl hexene-1. Little work has 
been carried out on the C, fraction. 

As a preliminary to a systematic investigation of the C, fraction, samples 
were obtained by fractionation in a Podbielniak column. The light 
hydrocarbon fraction from which it was derived came from the treatment 
of blue water-gas on a cobalt catalyst. A C, sample was treated with 63 
per cent. sulphuric acid according to the method of Dobryanskii.5> The 
absorption, assumed due to isobutylene, was 14 and 16 per cent. in two 
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cases. According to the method of Frey and Yant,®* in which absorption 
in a mixture of 2 vols. H,SO, (sp. gr. 1-84) and 1 vol. water is measured, 
16 and 19 per cent. isobutylene was found. 

The possibility of formation of aromatic hydrocarbons during synthesis 
may be considered to require an orientation of molecules at the surface 
of the catalyst which js ruled out by the conditions prevailing, particularly 
when dual attachment of the carbon monoxide molecule is involved in 
the mechanism of synthesis. No such special orientation is required for 
the production of branched-chain hydrocarbons of the paraffin and olefine 
series. 

In view of the importance of isobutane and isobutene as material for 
the production of motor spirit of high octane number and as material for 
chemical transformations for other purposes, a systematic investigation 
is planned to determine whether by modification of catalyst or in mode 
of operation the yield of these substances can be increased. 
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FRACTIONAL DISTILLATION OF BINARY MIX- 
TURES. SIMPLIFIED COMPUTATION OF THEO- 
RETICAL PLATES AND TRANSFER UNITS. 


By A. J. V. UNDERWoop, D.Se., M.I.Chem.E., A.M.I.Mech.E., 
F.R.LC., F.Inst.F., F.Inst.Pet. 


SuMMARY. 


Equations, in convenient form for numerical computation, are derived for 
the number of theoretical plates or transfer units required for the separation 
of binary mixtures of normal volatility. When the top or bottom product 
contains only a small amount of one of the components, computation can be 
simplified by using certain approximate equations which give satisfactory 
accuracy. By a simple transformation a calculation for a stripping column 
can be converted into one for a rectifying column. 


Ix a previous paper ' equations were derived for the number of transfer 
units required for the fractionation of a binary mixture, in which the 
components follow Raoult’s law and the relative volatility is constant. It 
was also shown there that these equations, as well as a modification of the 
equation published by Smoker ? for the number of theoretical plates, can 
be used to convert any separation with finite reflux into an equivalent 
separation with total reflux. Computation of a numerical example is 
however rather laborious when using these equations either for transfer 
units or theoretical plates. This is particularly the case for separations 
in which the top or bottom product contains only a small amount of one of 
the components. It will be shown that equations can be derived which 
offer considerable advantage in these respects. = 


NuMBER OF THEORETICAL PLATES 


By making use of a property of the rectangular hyperbole, an equation 
for the number of theoretical plates can be derived quite simply. In Fig. 1 


is shown part of a rectangular hyperbola, y = Boy which passes through 


the origin 0 and has asymptotes x = — B and y = A parallel to the axes 
of co-ordinates. A line through the origin, y = mz, intersects the hyper- 
bola again at Z. A line parallel to OX intersects the hyperbola at P and 
the line OL at Q. The co-ordinates of P are (x,, y,) and those of Q are 


Then - - - - = Gi) 


If the abscissa of L is l, then 
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Dividing equation (1) by equation (2), 
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This ratio is therefore constant for any pair of points such as P and Q 
determined by any horizontal line such as NM. From equation (2), 


B+l_ A 
B 
The slope of the tangent at a point on the hyperbola is given by 
dy AB dy A 
when x= 0. The constant _ is therefore equal to —S, where ¢ is 


the angle between the tangent to the hyperbola at 0 and the z axis and 6 
is the angle between the line OL and the z axis. 
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This is also evident from the following derivation : 


— 2) 
% 
As the line NM approaches the z axis and the points P and Q approach O, 


z, and 2» become small and — 1 becomes unity in the limit. 4% is 


0 1 
the tangent of the angle between the line OP and the z axis. When P 
and O coincide, the line OP becomes the tangent to the hyperbola at O. 


z ' is the tangent of the angle between OL and OX. As the ratio 
0 


is constant for any two points such as P and Q, it is equal to the value 
when P and Q coincide with O—i.e., 
Equation (3) may therefore be written : 


(4) 


ax 
Fig. 2 shows the usual equilibrium curve y = iti i 
a rectangular hyperbola with asymptotes parallel to the axes) and an 
operating line y = mz + 6 cutting the equilibrium curve at O’ and L, the 
abscisse of which are k, and k, respectively. P, and P are points repre- 
senting the compositions on two adjacent plates with abscisse z, and z, 
respectively. All co-ordinates now refer to the usual origin 0. 
It has been shown that 


(which is 


NQ 
QM tan¢d 
NP © tan 6 
PM 
— ky 
tang 
z,—k, tané 
k, — 2, 


¢ is the angle between the tangent to the equilibrium curve at O’ and the 
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and for = 
tan ¢ = 


where ¢ = 1 + (a — 1)k,. 
Also tan 6 = m. 
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Therefore, by successive multiplication : 


— 
(ra) 


which gives the number of theoretical plates required to pass from 2» to 2,. 
For total reflux, k, = 0, k, = 1,¢ = 1 and m = | and equation (5) reduces 
to the well-known form. n becomes infinite when z, = k, or xz = k,, 
representing the intersection of the operating line with the equilibrium 
curve for a rectifying column or a stripping column respectively. 

In the previous paper ! the equation derived by Smoker ? was converted 
into a form which may be written : 


a\ 2z,(1 — 2,) 
where = and zx, = 


Substituting these values in equation (6) then gives equation (5), which 
may also be written : 


(xy — — Zn) 


This is a convenient form for computation. k, and k, must first be found 
either graphically by drawing the operating line and noting its intersections 
with the equilibrium curve, or analytically by solving for the abscisse of 
the points of intersection from the equation : 


ak 
or m(a — . . (8) 


which is the equation originally given by Smoker.* 

It should be noted that k, is the lesser root of equation (8) and k, is the 
greater root. An operating line for a rectifying column will give a value 
of k, between 0 and 1, while &, will be greater than 1. For a stripping 
column k, will be negative'and &, will lie between 0 and 1. It should be 
noted that this convention is different from that used by Smoker, who 
does not use the root k,, and denotes as k that root of equation (8) which 
lies between 0 and 1 for both a rectifying and a stripping column. 

There is another simple proof of the relation given by equation (4a), 
which can be written in the form : 
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where d is a constant. The relation between x, and z,, as derived from 
the equilibrium curve and the operating line is : 


ar 


mm 


The constants k,, k, and d can be chosen so that the two equations are 
identical. The two equations can be written : 


— 1) + 2, (k, — dk.) + — dk,) + — 1) = 0 


and — 1) + 2,{b(a — 1) — a} + mz, + b= 0. 
The two are identical if 
d—1l &k,—dk, _k,—dk,_ k,k,(d —1) 
mia—1) 


These relations give three equations for determining k,, k,,andd. From the 
first, second, and third relations it is seen that : 

_ 

m(a—1) m+ (a —1)— 
m-+b(a—1)—a 


so that — (k, + &,) = 
b 
Also kk, = m(a — 1)" 
Consequently k, and k, are the roots of equation (8). 
Also (d — 1) = (a — 1)(k, — dk,). 

“T+ (@— 

Now 


{1s (a 1)k,}{1 (« 1)k,} =] (a 1)(k, + + (a 
+ (a — 1) — a} (a — 1) 


+ (a — me?” 


Although equation (5) or (7) is simple in form, its application to separations 
where the top or bottom product contains only a small amount of one of 
the components necessitates a high degree of accuracy in computation. 
Usually, for a rectifying column, z, = zp and z, = zp, and for a stripping 
column, x, = Zp and z, = Zy. When the top product contains only a 
small amount of the heavier component, the point L is close to the point 
C in Fig. 2 and k, is very nearly equal to 1. Also z, is very nearly equal 
to 1, so that the term (k, — zy) in equation (7) is the difference of two 
nearly equal quantities, which must therefore be computed with a high 
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degree of accuracy to avoid excessive error in the difference. Similarly, for 
a stripping column producing a bottom product containing a small amount 
of the lighter component, the point O’ will be close to O. &, will be small 
(and negative), and z,, will also be small, so that the term (x, — k,) will be 
small. Also for reflux ratios close to the minimum reflux, it will be clear 
that, for a rectifying column, the term (2, — k,) may be small, while for a 
stripping column the term (k, — 29) may be small. Exactly the same 
practical difficulty in computation arises if use is made of equation (6) 
or of the equation originally given by Smoker 2, namely : 


z,{1 — 1) 


im a — me* 


where 2’ = z — ky. 

The aioulty also arises when a graphical method such as that of McCabe 
and Thiele * is used, and necessitates drawing on a large scale those parts 
of the diagram where the points representing plates are close together. 
The nomograph published by Bisesi ‘ will also give insufficiently accurate 
results unless constructed on a large scale. 

For example, when a mixture of two components with relative volatility 
2-5 is fractionated with a reflux ratio of 3 to give a top product of composi- 
tion zp = 0-99, the value of k, is 1-0071 and the difference k, — zp = 
k, — z, = 90-0171. For three significant figures in the difference, five 
significant figures are required for the value of k,. 

The particular cases of a top product containing only a small amount of 
heavy component and a bottom product containing only a small amount 
of light component are of frequent occurrence in practical problems. For 
these cases simplified approximate equations can be derived which give | 
results of satisfactory accuracy with less laborious computation than the 
exact equations (7) or (9). Thomson and Beatty © have made a careful 
study of various methods of calculating the number of theoretical plates for 
such conditions. They have shown that, by assuming the equilibrium 
curve to be linear near the ends of the z, y diagram, reasonably accurate 
results.are obtained so long as the calculation is limited to moderate 
concentration ranges in these regions. For other cases, where there is a 
wide range between the feed concentration and the top or bottom product 
concentration, they do not recommend the use of that method. The 
simplified equations derived in the following are applicable to these 
conditions. 


NuMBER oF THEORETICAL CoLUMN. 
Put z) = tp = 1 — 8, where 3 is a small quantity. . 
R 
Substituting this value of b in equation (8) gives : 


m(x — 1)k® + {2m — ma — 1 — (a — 1)(1 — m)d}k + 
(1 — m)(1 — 8) =0 
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Solving this equation, neglecting powers of 8 higher than the first, gives : 


a(l — m)s 
and 
1 m(a—1)” m(ma— ma —1)l” (12) 


which gives the small term (k, — 2) in equation (7) with satisfactory 
accuracy by a slide-rule calculation. 

Equation (10) can also be readily solved by noting that k = 1 is an 
approximate solution, since it is satisfied by k = 1 when3 = 0. A second 
approximation to a solution is obtained by putting k = 1 + yf, where y is 
small, and neglecting small quantities of the second order, which gives : 


— m)s 


pat — SF an and k, = 1 + — a as before. 
Also, since k, and k, are the roots of equation (10), 
(1 — m)(1l — 3) 
—1) 
l—m («a — 1)8 


Of the four terms on the right-hand side of equation (7), (k, — 2) 
or (ky — xp) is given by equation (13). Also z, — k, = rp — k, (where 
xp is the feed composition) and is obtained from equation (12). If the 
feed is not composed of liquid at its boiling point, x» should really be 
replaced by 
(R + + (¢ — 
R+q 4 


where q is the heat required to vaporize 1 mol of feed divided by the molal 
heat of vaporization. In many cases this will be an unnecessary refine- 
ment unless g is quite small. 

The term (xz) — &,) is equal to (xp — k,) and ng term (k, — z,) is equal 
to (k, — 2p). 

It may be noted that 


= 


ma—1 — 1)8 
m(x—1) m(ma — 1) (14) 


In many cases no substantial error will be involved if the second term 
containing 3-is neglected. 


Also k, — tp = 1 — 2p + 


a(l — m)é 
ma — 


(15) 


and in many cases the term involving 8 can be neglected. 


Since 


U 

| An 

Here 

be for 

comp 

neces 

Fre 

subst 

log 3 

A 

meth 

comp 

ment 

Ale 

$ in 

equat 

This 

cent., 

error 

calcul 

In 

Since 

Equa' 

m( 

When 

small. 


the 


lal 


UNDERWOOD : FRACTIONAL DISTILLATION OF BINARY MIXTURES. 233 


An approximate value of — can be obtained to facilitate calculation. 
2(1 — m)(a — (16) 


ma —1 


me? ~ m{l + (a — 


Here again the term involving 3 will be negligible in many cases. It will 
be found that the use of equations (11) to (16) greatly facilitates numerical 
computation and avoids the use of five-figure logarithms which are generally 
necessary when using equations (6), (7), or (9). 

From the value of (k, — z,) given by equation (13), it will be seen by 
substituting it in equation (7) that is approximately a linear function of 
log 8 when m and « are constant. 

A comparison of the results given by the exact and the approximate 


ma{1 + 


methods is shown in Table I for various reflux ratios and top-product 


compositions, with « = 2-5 and rz» = 0-4. . It will be seen that the agree- 
ment is very satisfactory. 

A less accurate equgtion can be obtained by neglecting the terms involving 
8 in equations (12), (14), (15), and (16), giving, by substitution in 


equation (7) : 
n . log ma = log |- ime 


This equation will give substantial errors, of the order of perhaps 20 per 
cent., avhen used for reflux ratios near the minimum, mainly due to the 
error in the term (xy — k,). It may, however, be useful for certain rough 


calculations. 


NuMBER OF THEORETICAL PLATES—STRIPPING COLUMN. 
In this case zy is small and will be denoted by e. 


_ RP+qF _  — Way 


ba b = (1 — m)zy = (1 — m)e. 
Equation (8) then becomes : 
m(a — 1)k2 + {m + (a — 1)(1 — m)e — + (1— mje =0 . (18) 


When ¢ = 0, one root of the equation is equal to zero, and therefore k, is 
small. Neglecting small terms of the second order in equation (18) gives : 
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A second approximation for k, is obtained by putting 
a — 


m 


where ¥ is small. Substituting in equation (18) and neglecting small 
terms of the second order gives : 


2" 1) m(a—m) m(a—1) (a — 
The values of k, and k, given by equations (19) and (21) can also be obtained 


by solving equation (18) in the usual way and neglecting small terms of 
the second order. 


The approximation for nat 


a 2(a — 1)(1 — m)e 
_ (a— Ibe 


Equations (19), (21), (22), and (23) are then used for computation of 
equation (7). 

A comparison of the exact and approximate methods is given in Table II. 
As before, « = 2-5 and ry = 0-4. The results are given for various values 
of S, where S is defined by the operating line for a stripping column : 

S+1 


(S + = Sy + zp and m= ——. 


For a stripping column on which a rectifying column is superimposed : 


W_ (R + — (R + g)tw — (1 — 
Zp — Tp 


As previously pointed out by the author,* the quantity S may be termed 
the “‘ revaporization ratio.” It represents the number of mols of vapour 
returned by the reboiler to the stripping column per mol of bottom product 
withdrawn. It is thus exactly analogous to the reflux ratio in a rectifying 
column which represents the number of mols of condensate returned by 
the condenser to the rectifying column per mol of top product withdrawn. 
As a more convenient name for the quantity S, “ reboil ratio ” is suggested. 

It will be seen from Table II that the agreement between the exact and 
the approximate methods is very satisfactory. 

When equation (7) is used for a rectifying column, putting z, = xp and 
Z_ = Xp, the value of n obtained is the number of plates in the rectifying 
column, including the feed-plate. For a rectifying column used alone, 
mounted on a &till pot, the pot would be included as a plate instead of the 
feed-plate in the value of n. When equation (7) is used for a stripping 
column, putting x) = x, and x, = 2, the value of n obtained is the number 
of plates in the stripping column below the feed-plate and includes the still 
pot counted as one plate, For the common case of a rectifying column 
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superimposed on a stripping column the total number of actual plates will 
be the sum of the two values of n less one, the still pot being excluded. 

When « is not quite constant, it will usually be desirable to follow the 
suggestion of Thomson and Beatty ® and use the value of « corresponding 
to the region of the equilibrium curve where the plates are closest together. 
This will usually be at one of the ends of the equilibrium curve. 

A less accurate equation for the stripping column can be derived by using 
the value of k, given by equation (20) and neglecting ¢ in equation (22) 
and noting that ; 

Ly — ky = tp — ky = xp approximately 
and ky — 2, = k, — ¢ = k, approximately. 


i om 


This equation is subject to the same limitations that were pointed out in 
connection with equation (17). 


Then n. log = 


NuMBER OF TRANSFER UNITS. 


It has previously been shown ! that the number of transfer units is given 
by the equation : 


x 1 — 2, | 
ky 
where = met and ont 
Therefore 
me 


The equations previously derived for k,, k,, ete., can be used for the com- 
putation of equation (26) exactly as for equation (7). 


Mrximum Reritvux Ratio AND Minimum NuMBER OF PLATES— 
CoLuMN. 


The minimum reflux ratio, Ro, is given by 


Ry = — 
Yr — 
which reduces to 


xp(1 — zp)(x — 1) 


As a first approximation, which is sufficient for most purposes 
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The minimum reflux operating line with slope m, = —s intersects the 
0 
equilibrium curve at x = x,, which is the corresponding value of k,. 


Then, from equation (12) : 


— 
et 


mala — 1) — 1) 
an 
_ l-m l—m 


— 1) m(« — 1) 
From equations (11) and (29) : 


_ — a(l — m) 1 — 


ma—1 (m*« —1)8 
m(a — 1) ~ m(ma — 1) 


ma — 1 ‘ 
approx. . . . . (32) 
Then, substituting from equations (30), (31), (32), and (13) in equation (5) 
gives : 
a \* (ma — 1)*(m ga — 1) 
For total reflux, the minimum number of plates, n,, is given by 


— 2p) 


or 


From equations (33) and (34) : 


_ me — mp 


As an approximation from equation (16), 
= ma. 


a" (m—m,)(« — 1)? 


This equation can be used for plotting a curve showing the variation of 
n with m or R. 
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Resor Ratio anp Mintmom NuMBER OF PLATES— 
Srrieprinc CoLumn. 


As pointed out earlier in this paper, the slope, mg, of the operating line 
for the minimum reboil ratio, So, is given by 


—So+!_ yr — 


So tp — ty 
and = 
Yr — 


Noting that z, is the value of k, when m = mp, the following approximate 
equations can be derived by the same procedure used for the rectifying 


column : 


zp —k, = 
by — (41) 


Substituting in equation (5) from these equations and equation (23) : 


(a — m,)(a — m)? 
For total reflux, 

a —Zy) 


me*/ — m)*(m, — 1) 

As an approximation from equation (22), a 
(5) 

Then (45) 


Rerivux Ratio anp Resor Ratio. 


It has been pointed out that there is a close analogy between reflux ratio 
for a rectifying column and reboil ratio for a stripping column. It can 
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further be shown that, by suitable transformation, the calculation of 

plates for a stripping column can be made by using methods of calculation 

relating to a rectifying column. Fig. 3 shows an operating line for a 

stripping column, which cuts the equilibrium curve at D and F and the 

diagonal OB at E (x = zy). Its inclination to the z axis is given by 
S+1 


m = tan 0 = —_— If Fig. 3 is rotated through 180° about the diagonal 


B 


Fic. 3. 


AC, the result is Fig. 4. The equilibrium curve is unchanged since it is 
symmetrical with respect to a line joining A and C. The stripping-column 
line DEF of Fig. 3 has become a rectifying-column line in Fig. 4. It now 


has an inclination of (§ 0) to the horizontal axis and its slope is 
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Taking B as the new origin in Fig. 4, the co-ordinates of the point Z are 
(1 — 2p). If 2’ and y’ are the new co-ordirfates in Fig. 4, the equation of 
the line DEF is 


(46) 


This now represents an operating line for a rectifying column with a reflux 
ratio S and a top product (1 — zw). The calculation of plates for a stripping 


C - 


A 


Fie. 4. 


column can therefore be carried out by using the same equations as for a 
rectifying column, using S instead of R and (1 — zy) instead of zp. It 
ghould be noted that if z, is the feed composition for the stripping column 
in Fig. 3, the corresponding feed composition for the rectifying column 
in Fig. 4 is (1 — yp) = 2’. 

The above results can also be demonstrated analytically. The relation 
between the (zx, y) system of co-ordinates of Fig. 3 and the (z’, y’) system of 
co-ordinates of Fig. 4 is ) 


z=l—y andy=1-—2’ (47) 
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These equations indicate that the mol fraction of the lighter component 
in the liquid at any point in Fig. 3 becomes the mol fraction of the heavier 
component in the vapour at the corresponding point in Fig. 4. Writing 
the equilibrium curve in the form 


l-y 
it is obviously identical in form when transformed by using equations (47). 
Similarly the operating line 
= Sy+2y 
becomes (S + 1)y’ = Sz’ + (1 — zy). 
For a rectifying column superimposed on a stripping column there is a 


simple relation between R and S when the top and bottom products contain 
only a small amount of the other component. 


RP +qF —W_ (R + — (R + Q)tw — (1 — 
W tp — <p 


Since S = 


and zy is small and zp is approximately equal to 1, 


1 — zp 


When the feed is at boiling-point, g = 1 and then 


1 — 


This is also readily derived by considering that the top product is approxi- 
mately the lighter component of the feed while the bottom product is the 
heavier component. For 1 mol of feed there are zy mols of top product 
and (1 — zy) mols of bottom product. There are (R + 1) zy mols of 
vapour passing up the column and the number of mols of vapour per mol 


of bottom product are 
(R + 


1 — zp 
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Rectiryine CoLtumN—SumMARY OF RESULTS. 
ing 
Zp = 0-99. Zp = 0-97. Zp = 0-95. 
Appr. Exact. Appr. Exact. Appr. Exact. 
' 
R=158 . 
7). e- . | 0-361 0-3614 0-343 0-3447 0-324 0-3291 
1-0147 1-0145 1-0441 1-0424 1-0735 1-0691 
‘ - | 16354 1-6353 1-6922 1-6898 1-7490 1-7431 
(zy —&,) . | 0-039 0-0386 0-057 0-0553 0-076 0-0709 
(ky 2p) | 0-0247 0-0245 0-0741 0-0724 0-1235 0-1191 
4 . | 12-20 12-23 8-70 8-75 6-85 6-89 
3a R= 
in . | 06-3258 0-3259 0-3108 0-3120 0-2958 0-2989 
1-0125 1-0124 1-0375 1-0364 1-0625 1-0595 
a . - | 16917 1-6916 1-7417 1-7401 1-7917 1-7875 
(zp —&,) . | 00742 0-0741 0-0892 0-0880 0-1042 0-1011 
(k; — Zp) . | 0-0225 0-0224 0-0675 0-0664 0-1125 0-1095 
. | 10-45 10-46 7-65 7-72 6-19 6-29 
R=3 
-~ . | 0-2184 0-2184 0-2108 0-2111 0-2032 0-2040 
1-0071 1-0071 1-0214 1-0211 1-0357 1-0349 
= 1-8912 1-9232 1-9226 1-9558 1-9543 
(zp —k,) . | 0-1816 0-1816 0-1892 0-1888 0-1968 0-1960 
(ky —2p) . | 00171 0-0171 0-0514 0-0511 0-0857 0-0849 
- 3 . | 7-87 7-87 5-93 5-94 4-98 5-00 
R=4. 
_ . | 0-1642 0-1642 0-1592 0-1593 0-1542 0-1545 
1-0050 1-00499 | 1-0150 1-01488 | 1-0250 1-02464 
ci- = ‘ . | 2-1020 2-0119 2-0360 2-0357 2-0600 2-0595 
he (zy —k,) . | 0-2358 0-2358 0-2408 0-2407 0-2458 0-2455 
ct (ky 2p) . | 0-0150 0-01499 | 0-0450 0-04488 | 0-075 0-07464 
of ss . | 7-08 7-08 5-39 5-39 4-56 4:57 
ol 
* 
Stripping CotumN—SumMMARY OF RESULTs. 
ty = 0-01. Zw = 0-03. tw = 0-05. 
Appr. Exact Appr. Exact. Appr. Exact. 
S = 2. 
. . |—0-005 —0-00491 |—0-015 —0-01423 |—0-025 —0-02296 
. | 0-4527 0-4527 0-4694 0-4687 0-4861 0-4841 
a ‘ - | 1-6917 1-6916 1-7417 1:7402 -| 1-7917 1-7840 
(ky — zp) . | 00527 0-05269 | 0-0694 0-06867 | 0-0861 0-08407 
(tw —k,) . | 0015 0-01491 | 0-045 0-04423 | 0-075 0-07296 
. 5 . | 10-31 10-32 7-33 7-39 5-75 5-87 
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Srrierine CoLuMN—SuMMaARY OF (contd.). 


ty = 0-01. tw = 0-03. tw= 0-05. 
Appr. Exact. Appr. Exact. Appr. Exact. 
S = 3. 
| —0-002857 |—0-002836 |—0-008571 |—0-008349 |—0-01429 |—0-01368 
0-5887 0-5887 0-5994 0-5992 0-6101 0-6095 


1-9554 1-9545 


1-9230 


1-9232 


1-8911 1-8910 


(ky — 2p) . 0-1887 0-1887 0-1994 0-1992 0-2101 0-2095 
(tw — &,) 0-01286 0-01284 0-03857 0-03835 0-06429 0-06368 
n 7-17 7-17 5-22 5-22 4-24 4-24 


S = 4. 
ky - |—0-002 —0-001991 |—0-006 —0-005897 |—0-010 —0-009714 
0-6707 0-6707 0-6787 0-6786 0-6867 0-6864 


2012) 20120 | 2036 2-0359 | 2-060 2-0595 
(ky —zp) .| 02707 | 0-2707 | 0-2787 | 0-2786 | 0-2867 | 0-2864 
(zw —&k,) .| 0-012 0-01199 | 0-036 0-03590 | 0-060 0-05971 
Oa’ 6-30 4-59 4-60 3-76 3-76 


OBITUARY. 


MR. HIRAM J. HALLE. 


HreaM J. Hate, President of Universal Oil Products Company, died on 
29th May at his home, Pound Ridge, New York. His health had been 
failing for some time. 

Mr. Halle was born in Cleveland, Ohio, seventy-seven years ago. He 
entered petroleum refining at the time the automobile was coming into 
mass production, and for more than twenty-five years was an outstanding 
figure in the growth and development of the refining industry. 

He became president of U.O.P. in 1915, and under his leadership, the 
Universal organization took an important part in the invention, develop- 

, ment and servicing of processes which have revolutionized the refining 
industry, starting with the Dubbs cracking process and progressing through 

continuous research and development work to the perfecting of a number 

of catalytic refining processes. 
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CORRIGENDUM. 
JuLY JOURNAL. 


“ The Correlation of Cetane Number with other Physical Proper- 
ties of Diesel Fuels.” 

P. 196. Table 1. Column 3. Delete “‘ Outside’ before + 5. 
P. 197. Delete “(7)” before equation. 
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